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Galileon gravity offers a robust gravitational theory for explaining cosmic acceleration, having a 
rich phenomenology of testable behaviors. We explore three classes of Galileon models - standard 
uncoupled, and linearly or derivatively coupled to matter - investigating the expansion history 
with particular attention to early time and late time attractors, as well as the linear perturbations. 
From the relativistic and nonrelativistic Poisson equations we calculate the generalizations of the 
gravitational strength (Newton's constant), deriving its early and late time behavior. By scanning 
through the parameters we derive distributions of the gravitational strength at various epochs and 
trace the paths of gravity in its evolution. Using ghost-free and stability criteria we restrict the 
allowed parameter space, finding in particular that the linear and derivative coupled models are 
severely constrained by classical instabilities in the early universe. 

I. INTRODUCTION 

The origin of the current cosmic acceleration is a fundamental mystery, even to whether it arises from a new field or 
a change in the laws of gravity. Scalar-tensor gravity is an intriguing possibility, partaking of elements of both classes 
of explanation. Those theories that do not have a scalar potential avoid some naturalness problems associated with 
this, and those which involve a symmetry or geometric origin moreover sidestep difficulties with high energy physics 
corrections. 

A class of theories possessing all these desiderata is Galileon gravity [H-U, involving a shift symmetry in the field, 
while the field itself can be viewed as a geometric object arising from higher dimensions and entering in invariant 
combinations that assure second order field equations and protect against various pathologies. The derivative self 
couplings of the field also screen deviations from general relativity in highgradient regions (e.g. small scales or high 
densities), thus satisfying solar system and early universe constraints 0-0] ■ The galileon field, and closely related 
models [a4l(| , have been extensively studied in both theoretical fl7l - l25| and observational (2fj| - [33j contexts (and also 
with respect to inflation; (33-l36j). 

Galileon gravity thus merits investigation in detail, from both the expansion properties (homogeneous background) 
and structure growth (inhomogeneous perturbations) perspectives. Regarding the expansion, we will particularly 
be interested in the early universe corrections to radiation and matter domination, the effective equation of state 
evolution, and late time attractor behavior, especially to a dc Sitter state. For the behavior of linearized perturbations 
(concentrating on subhorizon scales), we explore the form of the modified gravitational strengths (Newton constants) 
in the generalized Poisson equations, looking for characteristic features in their evolution. 

The parameter space of the theory will be constrained by general physical considerations, such as positivity of the 
effective energy density and freedom from ghosts (Hamiltonian bounded from below). Observationally, we ask that 
the radiation and matter eras remain intact and the effective dark energy density is of order 3/4 the critical density 
today. These conditions still permit considerable diversity in the gravitational behavior, which is quite interesting to 
explore. 

In Section [TT] we briefly review the standard Galileon model and include terms with linear and derivative coupling 
to matter, setting up the field equations and modified Friedmann equations for a homogeneous, isotropic universe. 
Section [IIII sets up the linear perturbation theory and the generalized Poisson equations, deriving expressions for the 
gravitational strength modifications. Solutions for the expansion history for the uncoupled, linear-, and derivative- 
coupled models are presented in Section IIV1 including discussion of early and late time attractors. Evolution of the 
gravitational strengths is studied in Section [V] along with the no-ghost and stability conditions. We summarize the 
results and conclude in Section IVTl 



II. CLASSES OF GALILEON GRAVITY 



Galileon gravity is a scalar field theory containing nonlinear derivative self couplings. The action for the scalar field 
7r is invariant under Galilean symmetries tt — > tt + c + b^x^ in the absence of gravity, where c and are a constant 
scalar and four vector respectively. The four dimensional action that preserves these symmetries contains five unique 
terms, consisting of scalar combinations of d^ir, d^d^ir and On. 
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These models share many similarities with the decoupling limit of DGP gravity [37H39| and exhibit Vainshtcin 
screening that restores the theory to the general relativity behavior under certain conditions. Galileon gravity has 
field equations that contain at most second derivatives of the scalar field, and hence can be free of ghost instabilities. 
In addition, since these second derivative terms appear linearly, the Cauchy problem can be well defined. 

When gravity is introduced one must necessarily break the Galilean d^ir — > 
covariantized formulation of the Galileon model has been constructed in 
symmetry tt — > tt + c, and the Galilean symmetry is softly broken. 

The covariant Galileon action can be written as 
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where ci^s are arbitrary dimcnsionless constants, g is the determinant of the metric, M v \ is the Planck mass, and R 
is the Ricci scalar. The Galileon Lagrangians are given by 



A = Af 3 7r, C 2 = (V^ttXW), £3 = (□7r)(V^)(W)/Af 3 

U = (V^)(W) [2(Ott) 2 - 2tt,^ iv - E(V„7r)(W)/2] /M 6 

£5 - (V M 7r)(V"7r) [(Dir) 3 - 3(Dtt)tt.^tt^ + 2V - 67^7^7^^] /M 9 



(2) 
(3) 
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where M is the mass dimension taken out of the couplings to make the c's dimcnsionless; without loss of generality 
we can take M 3 = M p iH$. 

The Galileon invariance severely restricts the form of the action, although further freedom remains in the coupling 
between tt and matter. For example, the weak field limit of higher dimensional brane models typically involve a 
coupling of the form 7rT, where T is the trace of the energy-momentum tensor; such a term explicitly enters into the 
effective scalar field action in the decoupling limit of the DGP model [43[ . We also consider a derivative coupling 
TVdpTvdvTv, which has other interesting properties and origins. We will set the tadpole term to zero, i.e. c\ = 0, since 
we are interested in the effect of derivative self couplings on the growth and expansion histories and wish to avoid an 
explicit cosmological constant. 

Thus, the action can be written as 
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For computational purposes we find it is easier to work in the Jordan frame, where the explicit coupling between 
7T and matter is removed via a metric redefinition. In Appendix [X] we perform these metric transformations in the 
weak field limit, absorb like terms by renormalizing the c's, and promote the final actions to their full, non-linear 
counterparts. Here we state the result: 
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Varying the action with respect to the fields g^u and n, we obtain the general field equations that are exhibited in 
Appendix [B] 

In what follows we analyse three classes of models: the standard uncoupled Galileon, where cq = cq = 0, the 
linearly coupled Galileon, where cq = 0, and the derivative coupled Galileon, where cq = 0. Note that derivative 
coupling resulting in such G^ u d^dy-K terms has been found of interest for dark energy in (3j4r| and inflation in 
|47j |. and has ties to higher dimensional, vector, and disformal gravity theories. 

To derive the cosmological evolution of these classes of Galileon gravity, we use the Robertson- Walker metric for a 
homogeneous, isotropic spacetime, 



ds 2 



-dt 2 



5ijdx l dx° , 



(7) 



where a is the cosmic scale factor and we assume spatial flatness for simplicity. For the homogeneous evolution 
7T = 7r(a) and there are two independent field equations. 

We can write the (i, i) Einstein equation and tt dynamical equation in terms of dimensionless variables y = ir/M p \, 
H = H / H , and x = tt' /M„\, where the Hubble parameter H = a/a and a prime denotes dj 'din a, as a set of dynamical 
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equations 
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ct = 2(1 - 2c y)i? - 2c i?a; + 2c 3 H 3 x 3 - 15c 4 H 5 x 4 + 21c 5 H 7 x 5 + 6c G H 3 x 2 
A = 3(1 — 2c y)H 2 - 2c H 2 x 

-2c 3 H 4 x 3 + + ^ + ^c 4 i/V - Oc^ 5 - c G i/V 

2 a 4 2 

w = -2c i? 2 + 2c 3 H 4 x 2 - 12c 4 H 6 x 3 + 15c 5 H s x 4 + Ac G H 4 x . 
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As a check of our numerical solutions we verify that the redundant Fricdmann equation is satisfied at all times 
during the evolution: 



(l-2c y) H 2 



^mO . ^rO 



C2 fr2 2 



15 



+ 2c a H z x + -fH'x" - 2c 3 iTV + —c 4 H°x* - 7c 5 ff 8 a: D - 3c G £TV , 



(17) 



where £l m o and ^r-o are the present matter and radiation energy densities, respectively, in units of the critical density. 
It is also useful to write the energy density and pressure of the scalar field: 



H 2 M 2 , 



6c Q H 2 x + ^H 2 x 2 - 6c 3 H 4 x 3 + yc 4 iJ 6 i 4 - 21c 5 H 8 x 5 - 9c G H 4 x 2 



H 2 M 2 vX 



= -co 



4H 2 x + 2H (Hx)' 



^H 2 x 2 + 2c 3 H 3 x 2 (Hx)' - c 4 



^H 6 x 4 + 12H 6 x 3 x' + 15H 5 x 4 H' 



+3c 5 H 7 x 4 (pHx 1 + IH'x + 2Hx) + c G [6H 3 x 2 H' + AH 4 xx' + 3H 4 x 2 ] 
from which one can define an effective equation of state parameter for the Galileon: w = P^j p^. 



(18) 



(19) 



III. INHOMOGENEOUS PERTURBATIONS AND MODIFIED POISSON EQUATIONS 

While the previous section set up the background dynamics for the cosmological evolution, we also need the 
inhomogeneous equations to investigate the effects of perturbations on the cosmic growth history and gravitational 
modifications to the generalized Poisson equations. We here consider linear perturbation theory for scalar modes in 
the subhorizon limit. In the Newtonian gauge the perturbed metric is 

ds 2 = -(1 + 2i')dt 2 + a 2 (l - 2(f))S lj dx l dx j . (20) 

The linearized Einstein and scalar field equations are exhibited in Appendix [Cj These lead to the modified Poisson 
equations 

VV = W g|»"- a „ (21) 

v V= ^gky (22) 

V 2 W + ct>)= 87raG «% Pm 6 m , (23) 
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where we describe the effect of the Galileon on subhorizon density perturbations through modifications G e g of the 
gravitational strength, or effective Newton's constant. These will be functions of the background quantities H, y, and 
their derivatives. 

Note that G^g is central to the growth of cosmic structure and is called V in [48j| and p in [U, while G% + * ] 
is important for light deflection and the integrated Sachs- Wolfe effect, and called Q and S in those two references, 
respectively. Translation tables between other parametrizations are given in [3 [Hlj . One can regard the ip Poisson 
equation as governing nonrelativistic geodesies and the ip + 4> Poisson equation as governing relativistic geodesies in 
the inhomogeneous spacetime [Hi], [HH . 

The modified gravitational strengths are given in terms of Newton's constant Gn by 

g<&> = „ ,„ ,. — — g n ( 24 ) 

G it ] = ZTTZTZ ' J'tL Z TTTT (25) 

G (W) = - uv -o , ^ -,v ^ Gn 

K5 (K4K1 — K5K3J — K4 (K4K6 — K5K1) 

where ki_6 are functions of the background variables: 

Ki = -6c 4 H 3 x 2 (h'x + Hx' + + 2c G (HH'x + H 2 x' + H 2 x) - 2c + c 5 H 5 x 3 (l2Hx' + IhH'x + 3Hx) 

k 2 = -y + 6c 3 H 2 x + 3c G H 2 - 27 Ci H 4 x 2 + 30c 5 H 6 x 3 

k 3 = -(1 - 2c y) - yifV + c G H 2 x 2 - 3c 5 H 5 x 4 (Hx' + H'x) 

k 4 = -2(1 - 2c y) + 3c 4 H 4 x 4 - 2c G H 2 x 2 - 6c 5 H 6 x 5 
k 5 = 2c 3 H 2 x 2 - \2c A H 4 x 3 + Ac G H 2 x - 2c + 15c 5 H 6 x 4 

k 6 = y - 2c 3 (H 2 x + HH'x + 2H 2 x) + c 4 (\2H 4 xx + l8H 3 x 2 H' + 13H 4 x 2 ) 
-cg (2HH' + 3H 2 ) - c 5 (l8H 6 x 2 x' + 30H 5 x 3 H' + 12H 6 x 3 ) . 

The general relativistic form of the Poisson equation is recovered if we set cq = 03-5 = c G = (02 can remain 
unspecified; in this limit we have a massless, minimally coupled scalar field which will not modify the effective Newton 

constants). In that case, G { ^> = G^ = G^ +4>) = G N . 



IV. EXPANSION HISTORY IN GALILEON COSMOLOGY 



We now evaluate how Galileon cosmology affects the cosmic expansion history relative to the standard ACDM 
paradigm, including early universe asymptotes, the onset of cosmic acceleration, and late time attractors, especially 
to a de Sitter state. 



A. Uncoupled Galileon: Co = = cq 



The first, simplest case corresponds to the absence of any direct coupling between the scalar field 7r and matter; 
the parameters are then 02,03,04,05. Such a model has been extensively studied in [30M32I] . Equations (j5])- (fTU)) 
describing the background evolution define an autonomous system, and it is straightforward to establish the existence 
of asymptotic fixed points. 

The only fixed points for matter and radiation domination correspond to (fi m , £l r , x) — (1, 0, 0) and (0, 1, 0) respec- 
tively. We therefore solve the 7r equation assuming matter and radiation domination to obtain the behaviour of p n at 
early times. For late times, we want cosmic acceleration and so we look for a late time de Sitter asymptotic state by 
finding the fixed points of the dynamical system (H,x), taking f2 m = O r = 0. 

The full numerical solutions are illustrated in Figure[T]for specific parameter choices, and agree asymptotically with 
the fixed point behaviors. We plot the time evolution of H 2 , p m /(3H 2 M 2 l ), p t /(3HqM^) and p n / '(3H 2 M 2 X ) in the 
left panel, and the effective dark energy equation of state parameter w in the right panel. 
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FIG. 1: [Left panel] The redshift evolution of H 2 (solid red), 87rGp m /3i?o (dotted blue), STrGpr/SHg (dot-dashed magenta) 
and 8ttGpk/3Hq (long dash black) for the uncoupled Galileon, with parameters C2 = —27.3, C3 = —12.8, C4 = —1.7, C5 = — 1. 
Here, the scalar field energy density slowly grows relative to matter at early times, as expected from the analytic results. The 
C5 term dominates the n dynamics up to z ~ 10, after which we observe a turnaround to approach the late time de Sitter point. 
[Right panel] The effective equation of state of the Galileon field. During radiation domination w — —1/4, in agreement with 
the analytic results; the model currently accelerates the expansion and approaches a late time de Sitter point. 



To calculate the behaviour of the tt field, at early times we look for solutions taking H ~ a 2 (during radiation 
domination) or H ~ a -3 / 2 (during matter domination). The field equation of motion (JSJ becomes 

_ -5c 3 H 4 x 2 + 30c 4 H 6 x 3 - 75c 5 H 8 x 4 /2 

(27) 



c 2 H 2 /6 - 2c 3 H 4 x + 9c 4 H 6 x 2 - 10c 5 ff 8 a; 3 

during radiation domination and 

, _ -c 2 H 2 x/12 - 7c 3 H 4 x 2 /2 + 45c 4 H 6 x 3 /2 - 115c 5 H 8 x 4 / 4 

X ~ ~ X + c 2 H 2 /6 - 2c 3 H 4 x + 9c 4 H 6 x 2 - Wc 5 H s x 3 ^ ' 

during matter domination. We note the existence of a solution x ~ H ~ 2 during both radiation and matter domination; 
this behaviour was first observed in (32| . However, solutions close to this fixed point at early times require an 
extraordinary fine tuning of initial conditions to ensure that the tt field energy density remains subdominant until late 
times, as the energy density of the tt field grows at a fast rate relative to matter and radiation (the tt field equation 
of state is w = —7/3 and w = —2 during radiation and matter domination respectively.) We are interested in the 
behaviour of the tt field for generic initial conditions. 

Since each C2-5 contribution to the scalar field energy density enters in the ratio of H 2 x relative to the previous 
one, whether the initial condition has X{ <C H^~ 2 or not determines the early time evolution. For example, if this 
holds then the lowest, c 2 term will dominate p n and x ~ a -1 , x ~ a -3 / 2 during radiation and matter domination 
respectively. Since the other terms then decay at a faster rate than the c 2 term, once c 2 dominates it will do so at all 
subsequent times. In this case the tt energy density will decay as a massless scalar field p„ ~ a -6 at all times (recall 
c 2 involves the canonical kinetic term), leading to no interesting dynamics. 

As we increase Xi the 03-5 terms will dominate. This is expected in the early Universe, where the model will be in 
the strongly (self)coupled regime. Note that the dynamics of p n is hierarchical: for example if the c 3 term dominates 
at a particular time then the 04,5 terms will remain subdominant at all subsequent times, however the c 2 term grows 
relative to the c 3 term and can eventually dominate. 

During the radiation epoch, we find the following behaviour of x and p^, depending on which of the C2-5 terms 
dominate: 

(29) 
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(32) 
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and during the matter epoch 
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The effective equation of state of the it field varies between w G [—1/4, 1/4] when the 03-5 terms dominate during 
these epochs. During radiation domination, the 7r energy density typically grows relative to the matter component, 
however this effect is generally small and does not require a significant tuning of Q„(zi) initially. For example, if 
the C5 term provides the main contribution to p^ during radiation domination then we have p v /p m ~ a 3 / 8 . This is 
the worst case scenario, and provides a limit on how large we can set p^ initially such that it remains subdominant 
at matter radiation equality (the exact constraint depends on our choice of Zi). Thus Galileon models possess the 
interesting property that they can belong to the class of early dark energy models, in the sense that their energy 
density is non-negligible at early times for generic initial conditions. (Indeed, the C4 dominating case maintains a 
constant fraction of the matter density during matter domination.) 

To calculate the behaviour of the n field at late times one can obtain the asymptotic dc Sitter fixed points using a 
dynamical systems approach based on Eqs. ([I])-®. The fixed points correspond to solutions of the algebraic relations 

A = —ujx (37) 
7 = -fix (38) 

that arc coupled polynomials for H 2 and x. For any values of the parameters C2-5 there is a double zero corresponding 
to Minkowski space; H = 0. In addition, one can show that there are at most three de Sitter, three anti de Sitter 
and six unphysical complex solutions for H . To see this, we can use the coupled parameter x = H x, in which case 
Eq. (|38|) reduces to the cubic polynomial 

15c 5 x 3 ~ I8C4X 2 + 6c 3X - c 2 = , (39) 

which has three (one) real solutions if its discriminant 

D = lUcjcl - 160c^c 5 + 360c 2 c 3 C4C5 - 288c 2 c| - l^c\c\ (40) 

is positive (negative). One can also state that for C5 < (as needed for positive energy density when eg dominates 
early) there is exactly one real solution for x whenever 5C3C5 — 6c 2 > 0. Note that this does not necessarily imply the 
existence of a de Sitter point, which can only be determined by solving Eq. (j3"7]) . In terms of x, the condition for the 
existence of a de Sitter vacuum state is 

3i? 4 = -6 C5X 5 + ^4X 4 -|x 2 >0. (41) 

Hence for any real \, there is at most one positive and one negative real H solution to Eq. (|4ip . and there are at 
most three distinct, real \ solutions to Eq. ([39| . The minimal number of de Sitter solutions is zero. 

Certain specific cases are worth mentioning. We require at least two non-zero c 2 _5 terms to admit a dc Sitter fixed 
point, and during the approach to the de Sitter asymptote, gcncrically all terms are of the same order. Restricting 
our analysis to the case that any two of the terms C2-5 are non-zero, we can succinctly write the conditions for 
positivity of p T at early times and the existence of a late time de Sitter point as: (02, C3) = (— , — ), (C2, C4) = (—,+), 
(c2, c.5) = (— , — ). In each case, there is exactly one de Sitter point (although in the C4 case the solution has multiplicity 
two). Note that the standard kinetic term must exist (or the no-ghost condition of Sec. IV Al is violated) and must 
take the non-canonical sign. 

Both the early time and late time analytic results are borne out by the numerical solutions of the dynamical 
equations. Throughout this article we take initial conditions at z% = 10 6 of Pn(zi) = 10 -5 p m (zi) and fix 0, m {z = 0) = 
0.24. The parameters, and behaviors, exhibited in Figure Q] are typical in the sense elucidated in Section IV Dl 



B. Linearly Coupled Galileon: co 7^ 0, cg — 

We now consider a linear coupling between the field 7r and the trace of the energy-momentum tensor (or Ricci 
scalar) (26l. l42l|. As mentioned, such a scenario is typical in higher dimensional braneworld models, and in particular 
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in the decoupling limit of the DGP model. We begin as before by examining the behaviour of the tt field during 
matter and radiation domination. We initially consider the DGP like case where C4.5 = and then extend to include 
the more general Galileon kinetic terms. 

In the decoupling limit, the DGP model can be written as a four dimensional effective action with a scalar field 
coupled to the Ricci scalar and containing kinetic terms of the form C2.3. With the parameters Co, 02,03 we expect 
to recover DGP-likc behaviour in the early universe, where p^ is subdominant and we can neglect the backreaction 
of the 7r field on the geometry. Indeed, if we choose initial conditions such that coy < 1, we find two approximate 
solutions during matter domination of the form H ~ a -3 / 2 , x = ±Ao-ff _1 , where Aq = c^/{Zc^). On this solution, 
Pt< ~ ±AcqAqH , the familiar modification to the Friedmann equation one would expect from a model mimicking the 
DGP model (here we choose the sign convention c , C3 > 0). 

On the "self accelerating branch" corresponding to x > 0, p„ > 0, one can also see that the tt field does not satisfy 
the no-ghost condition given by (see Sec. IV A[) 



(c 2 - 12c 3 H 2 x) (1 - 2c y) + 6 (c 3 H 2 x 2 - c f > . (42) 

At early times the first term on the left hand side dominates, and we can neglect the C2 and coy contributions, leaving 
just — 12c 3 H 2 x > 0. For positive C3, the sign of x dictates whether the field is a ghost; on the self accelerating branch 
the inequality is not satisfied and tt is a ghost as expected [HI, [H[ . 

We now introduce the 04,5 terms, and consider how the evolution of tt is modified. To exhibit the behaviour of x at 
early times, we write H ~ a~ 2 and H ~ a~ 3 / 2 for radiation and matter domination. But it should be observed that 
the standard behaviour of H is no longer guaranteed: from the Friedmann equation (|17p it is clear that to obtain the 
standard behaviours H ~ a~ 2 and H ~ a~ 3 / 2 we must impose Coy < 1 to avoid modification of the left hand side. 
We therefore set y — initially, at some redshift deep in radiation domination, and check that y remains small. 

During radiation domination, we find the same solutions ([2"9| - (l32l as in the uncoupled case due to the fact that the 
Co coupling term is subdominant in the x equation, so 

_ , -5c 3 H 4 x 2 + 30c 4 H 6 x 3 - 75c 5 H s x 4 /2 



c 2 H 2 /6 - 2c 3 H 4 x + 9c4,H 6 x 2 - 10c 5 H s x 3 
However during matter domination we find modified behaviour, 

, -c 2 H 2 x/12 - 7c 3 H 4 x 2 /2 + 45c 4 H 6 x 3 /2 - c H 2 /2 - 115c 5 H 8 x 4 /4 



c 2 H 2 /6 - 2c 3 H i x + 9c 4 H e x 2 - l0c 5 H s 



(44) 



The behaviour of x and hence p v depends upon the initial conditions imposed, and general analytic solutions do 
not exist. However for certain limiting cases one can construct explicit solutions. One example is when the c 3 term 
dominates, corresponding to the DGP like solution discussed above with p^j (H 2 M 2 ^) = AcqAqH. Analogous solutions 
exist for the cases where other terms dominate; note that again the contributions to p n are hierarchical - if the c n 
term dominates early, it can only be superceded later in matter domination by a c m term with < m < n. The 
matter dominated era solutions, depending on which of the C2-5 terms dominate alongside cq, are 

c 5 ,c : x = A Q H-V 2 p^^coAoH'/ 2 A = ( w = ~ (45) 



5 u u u Vl5c 5 

1 /3 

AoB-W P7! = 7 -c A H 2 / 3 A =(-^-) w = ~ (46) 



2 u u V 9c 4/ 21 



x 1/2 

CO 



c 4 , c : x 

c 3 ,c : x^AqH" 1 Pt ! =Ac A H A c ,= [j^J w = -1 (47) 

c 2 , c : x = A pn = 5c Q A H 2 A = — — w = -% , (48) 

c 2 5 

where we have chosen sign conventions such that x > 0. Note that because of the explicit coupling to the Ricci scalar, 
the standard continuity equation and p n ~ a - 3 ( 1 + tu ) do not hold. We also note that the last case is only viable if 
I c§ / C2 1 <C 1, otherwise we will pick up non- negligible corrections to the standard matter era due to the presence of 
the coy term on the left hand side of the Friedmann equation. We will show in Appendix [D] that in each of the other 
cases, the it field is a ghost whenever p^ > (this result does not depend on our choice x > 0). Therefore the linear 
coupled model is only viable when either the C2 term dominates the c m >2 terms at early times and |cq/c2 | 1, or 
when the c 2 _5 contributions to the it energy density are much larger than that from c . In the latter scenario the 
model will evolve approximately as in the uncoupled case during matter and radiation domination. 




FIG. 2: [Left panel] The redshift evolution of viable linear coupled Galileons is restricted to be near that of the equivalent 
uncoupled Galileon, except at late times. The late time behaviour of H 2 for the linearly coupled model (dotted black) is 
compared to the uncoupled model H 2 (solid red) with the parameter choices (02, C3, C4, C5, Co) = ( — 5.9, — 2, —1, —1,0.05) and 
(—6.0, —2, —1, —1, 0) respectively (C2 is adjusted to ensure Q n ,o = 0.76 in both cases). Note the Hubble parameter in the linear 
coupling case no longer asymptotes to a constant value. [Right panel] The effective equation of state for the linearly coupled 
7r field as a function of redshift. At early times w resembles the uncoupled case, e.g. Fig. [T] but deviates at later times and in 
the future has w close to but not exactly —1. 



At late times we no longer obtain a de Sitter fixed point; this is due to the fact that the field y now enters the 
Fricdmann equation explicitly. The de Sitter solution in the uncoupled case corresponds to H , x approaching constant 
values asymptotically, and this is no longer possible since then coy would grow. However, since y grows as In a for 
constant x, a near de Sitter state is possible, where w ~ — 1 for z < —0.5. From our numerical studies we find 
that typical deviations from w = —1 are less than ~ 5% for model parameters in the range C2-5 ~ (—10, 10) and 
c ~ (—1, 1) at z ~ —0.95 (though larger deviations are possible). 

The numerical evolution of the linearly coupled model is exhibited in Figure [5] for generic parameter values cq = 0.1 
C3 = —1, C4 = 1, C5 = —1, which fixes C2 = —10.0 so that SX % {z = 0) = 0.76. The high redshift evolution is necessarily 
similar to the uncoupled case, as just discussed, but the future evolution differs. We note that H and x continue 
to evolve in the future and do not approach constant values as in the uncoupled case, so there is no de Sitter state; 
nevertheless the equation of state w remains very close to —1. 



C. Derivative Coupled Galileon: cg 7^ 0, co = 



In this subsection we include the derivative coupling term with cq, and switch off the linear coupling term with 
Co- Recall that the kinetic coupling to the Einstein tensor, G^d^ird^iT, also arises in some higher dimension gravity 
theories and disformal field theories and by itself has interesting behavior involving cosmic acceleration [3, |45| . 

At early times, the cq contribution will dominate over the C2 one, and for reasonable initial conditions (i.e. that 
the energy density p w <C p m initially) x will satisfy i« 1 and hence the cg term will dominate over the C3 term as 
well. However the importance of the cq term relative to c^s depends on the specific initial conditions. If cg does not 
dominate, then we have the previous uncoupled case behaviors. If cg does dominate then during radiation domination 

cq : x ~ a 3 p n ~ cqcl~ 2 w = — 1/3 (49) 

and during matter domination 

cq : x — a 3/2 pn — c G a~ 3 w = . (50) 

Note that p^ grows relative to the matter energy density during the radiation epoch and scales with it for z < 10 3 , 
and therefore we must be careful to choose initial conditions such that p n remains subdominant until z ~ 1. In 
particular, we impose that < 2 x 10~ 2 fi m at matter/radiation equality. We exhibit the full behaviour in Figure^ 
numerically evolving the equations for the case 04 = 05= 0, C3 = —1 and cq = — 1. The transition from w = —1/3 
during radiation domination to w = for z < 1000 is clear. If we switch on the 04,5 terms then since those terms 
grow relative to the cg term during radiation domination the limiting cases of Eqs. (|49[) and (1501) will no longer be 
clearly distinguishable. 
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FIG. 3: [Left panel] The evolution of H 2 (solid red), 8nGp m /3Ho (dotted blue), 8irGp r /3Hq (dot-dashed magenta) and 
8ttGpt,/3Hq (long dash black) for the derivatively coupled Galileon, with C2 = —9.1, C3 = —1, C4 = eg = 0. [Right panel] The 
equation of state of the -k field as a function of redshift. The numerical evolution exhibits the analytic attractor results. 

At late times the derivative coupled Galileon goes to the de Sitter attractor. Note that in the model of [3] using 
derivative coupling, the c 2 term was forced to the canonical form, and so no stable de Sitter attractor was allowed, 
however here there is sufficient freedom to enable an asymptotic de Sitter state. 

V. EVOLUTION OF PERTURBATIONS AND GRAVITATIONAL STRENGTH 

Having established the behaviour of the homogeneous cosmological expansion for the classes of Galileon models, 
we now examine the behaviour of linearized perturbations. Since the derivatives of the 7r field couple to the metric 
potentials, we expect the field to modify the growth of matter density perturbations by introducing time dependent 
couplings between Sir, 4> and V'- The equations governing the behaviour of subhorizon perturbations are given in 
Appendix [C] As discussed in Section Mil the effects on growth of structure (and light deflection) are concisely en- 
capsulated in the gravitational strength modification functions G^g and G^ + ^ . It has been shown in Sl| that 
stringent constraints can be placed on the Galileon class of models by considering the time evolution of Gn- 

The "paths of gravity" describe the evolution of these functions, from the initial deviation from the general rela- 
tivistic value of Gn at high redshift to present day signatures and the eventual asymptotic values. It is important 
to stress the domain of validity of our equations; they are applicable to linear perturbations on subhorizon scales, 
assuming that dark matter and baryons can be treated as a common fluid. We consider the evolution of perturbations 
over the redshift range z < 500. We set initial conditions at z — 10 6 as = 10~ 5 p m , to ensure that at z ~ 500 the 
perturbations remain small and G e s is sufficiently close to its general rclativistic behaviour. At low redshifts, linear 
perturbation theory is only applicable for large scale modes. To model non-linear perturbations at the present time we 
would be forced to resort to simulations, which is beyond the scope of the present work. Significant deviations from 
the linearized behaviour will be expected in the non-linear regime, due to the presence of a cosmological Vainshtein 
screening effect. 

The perturbation equations for the tr field itself must also satisfy certain physical conditions for a sound theory, 
such as the absence of ghosts modes and Laplace instabilities. These will constrain the allowed parameter space, as 
can conditions placed on the sign and magnitude of G e ff. Before discussing the gravitational evolution we briefly 
discuss the no-ghost condition, later treating its details and the Laplace instability in Appendices [D] and |EJ 

A. No-Ghost Condition 

The equation of motion for the field perturbations takes the form of a wave equation, with the sign of the second 
time derivative Sy needing to be non-negative to avoid ghosts (negative kinetic terms causing the Hamiltonian to be 
unbounded from below). For the standard, simple case of a minimally coupled, canonical scalar field the condition is 
just positivity of the kinetic energy, i.e. c-i > when all other c's are zero. In the presence of nonlinear kinetic terms 
in the action, and couplings, the situation is more complicated and we present our analysis in Appendix [D] Here we 
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summarize the result: to be free of ghosts the theory must satisfy 

c 2 - 9 , - fi , (3c 3 H 2 x 2 + 6c G H 2 x-l8c 4 H i x 3 + fc 5 H 6 x' l -3c ) 2 . 

+ 6^ + ^-27^+30^x3 + 2 ^ _ 6(1 - 2coy) 6cgR2x2 +Qc4g 4 4 1 lgc6g6 / < 0. (51) 



B. The Thawing of Gravity at Early Times 

At early times, with the Galileon contributions to the energy density small, the gravitational strength functions go 
to the general relativistic values of unity. That is, the early universe behaves as in standard gravity and cosmology. 
As the Galileon energy density increases, the modifications to gravity grow; we can say that the theory thaws away 
from general relativity, in analogy to the thawing class of dark energy dynamics (where the scalar field moves away 
from a frozen, cosmological constant state). 

For the initial stages of thawing we can calculate both G e ff and the no-ghost condition analytically for each of the 
cases where a given term of C2_5,g dominates the n energy density. For the no-ghost condition, we find 

C2-- -fk<° ( 55 ) 

where p^ = p„/{HqM^ is the dimensionless Galileon energy density. In all cases, the condition p n > is sufficient 
to ensure that the field is not a ghost. Note that given the initial conditions such that p^ <C H 2 during matter and 
radiation domination, the first terms in each will always dominate. 
For the gravitational strength, we find 

(57) 
(58) 
(59) 
(60) 
(61) 



C5 : 


Gn 


= 14 


4a 

7 


c 4 : 


r (4>) 

U cff 

Gn 


= 14 




C3 : 


Gn 


= 14 


5 


C2 : 


r {<t>) 
Gn 


= 1 




cg : 


Gn 


= 1 - 





These analytic solutions are only valid during matter domination. Figure |4] exhibits the evolution of the deviations 
from general relativity G e ff/Giv — 1 for each of the effective gravitational strengths, from an early time where they 
share similar redshift evolution in the deviation from general relativity, to a shared late time de Sitter asymptote as 
discussed in the following section. 
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FIG. 4: The thawing of gravity away from its general relativistic strength, G^t /G 



1 (solid red), G^/Gn 



1 (long dash 

black), and G^ + *'/Gn — 1 (short dash blue) for model parameters (02, C3, C4, C5) = ( — 15.3, —5.73, —1.2, —1). All three functions 
share a common redshift dependence at early times and a common future asymptote, as discussed in the text. 



C. Gravity in the de Sitter Limit 

The other limit of importance is the late time approach to the de Sitter attractor that is present in the uncoupled 
and derivatively coupled models. In the de Sitter limit, we have H' , x' — > 0, and it is convenient to define 

(62) 
(63) 
(64) 
(65) 
(66) 

Note that generically all these terms are of the same order (unless a coefficient is zero) so we cannot study the limit 
in terms of a single dominant component, making the analysis more difficult but the behavior richer. 

In general (with one type of coupling) we have five parameters and two independent equations of motion. The 
background field equations then reduce to algebraic constraints: 



c = 


c 2 x 2 


D = 


c 3 H 2 x 3 


E = 


Ci H 4 x 4 


F = 


c 5 H 6 x 5 


A = 


c G H 2 x 2 



E=- (-2A + 8D - 3G - 10) 

f=1(d-a-°- 

3 V 2 



(67) 
(68) 



Eliminating E, F from G^g , we arrive at the de Sitter limit 

G e ff,dS : 



14 + 7 A + 3C - AD 



G 



N • 



(69) 



Interestingly, in the de Sitter limit = G^g = - a special property of the Galileon models that in this 



limit they have no gravitational slip - and we write them simply as G ff,ds- 

It is also useful to write the no-ghost condition at the de Sitter state (although the no-ghost condition must hold 
at all times) in terms of A,C,D. We find 



^G + 2£>-5A-10)^ 



-C + 2D- 11A - 10 + 2-^-2- 
2 8 + 3G - AD + 4A 



< . 



(70) 
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For the uncoupled case A — 0, the no-ghost condition corresponds to 

2+^<D<5+^ . (71) 
4 4 

Figure [5] illustrates the constraints imposed on the Galilcon parameter space by several conditions for physical 
viability, including the no-ghost condition, Laplace stability condition on the sound speed, and positivity of energy 
density p v , and their relation to the asymptotic dc Sitter value of the gravitational strength G ff,ds- For the uncoupled 
case note that C5 < 0, which is generally required in the early universe barring fine tuning, requires that both C3 and 
C4 must also be negative. 




FIG. 5: [Left panel] Dark red shaded region shows the area violating the no-ghost condition in the C-D plane, fixing C5 = 0. 
The straight lines show the cuts through the plane corresponding to special values c n — as labeled. The light cyan shaded 
region shows the Laplace instability c? a < 0, and the thin "bubbly" red region shows the superluminal condition c? a > 1 (see 
Appendix [E]). G > is also ruled out. [Right panel] For the cg = case the no-ghost region lies between the solid black lines. 
These boundaries correspond to the limits G e ff,ds/Gjv = ±1/2; the blue dashed lines give G e ff,ds/Gjv = ±1, and the red dotted 
lines give G e fr,ds/Gjv = ±5. Positive values are for the upper line of each set. G e ff,dS passes through infinity for 4D — 3G = 14. 
Note G e fF,ds/Gjv cannot lie between —1/2 and 1/2. Further constraints include positivity of scalar field energy density at high 
redshift, limiting the parameter space to the left of the long dashed, diagonal cyan line when the C5 term dominates, and the 
Laplace inequality (|E3I) . with shaded magenta regions obeying < < 1. 



Since after the equations of motion are applied one has G e ff as a function of three free parameters, it is difficult to 
make general analytic statements. We begin by considering restricted cases, and then carry out a more general scan 
through the full parameter space. 

First suppose that only one of the parameters C2-5.G is nonzero. We find that in no case is there a consistent 
de Sitter solution. Solutions are possible when, e.g., two or three C2-5,g terms are nonzero. For a more physically 
motivated restriction, one can look at parameter choices for which gravity returns asymptotically to general relativity, 
at least in the sense of G e ff,ds = 1- This requires the condition 

17 + 7A + 3G-4L> = 0, (72) 

so now our 5 dimensional parameter space is restricted to 2 dimensions and we can plot allowed regions, imposing as 
well the viability conditions. The no-ghost condition becomes 

3G-4^<- 39 - 7 ^ or 4>3G-4^>- 39+ 6 7 ^, (73) 

(although the second region violates Laplace stability, having (? s < 0) and is exhibited in the left panel of Figure |U 
Setting some c n = collapses the 2 dimensional space to a 1 dimensional line. 
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FIG. 6: [Left panel] Region of the Galileon parameter space that gives a de Sitter limit of G b s/Gn = 1 and has no ghosts 
or Laplace instability is in the unshaded region. The 2 dimensional viable region becomes a 1 dimensional line if further 
restrictions are imposed. Note the uncoupled case cg = is here always ghost free, though it is Laplace unstable, <? a < 0, for 
insufficiently negative C. [Right panel] The evolution of the effective gravitational strengths G^ (solid red), G^g (long dash 
black) and G^^ (short dash blue) for an uncoupled Galileon case restoring to general relativity in the de Sitter limit, with 
and (c2, C3, ca, eg) = ( — 12.688, —6.239, —2.159, —1). Note the strengths can be significantly different from each other. 



D. The Paths of Gravity 



Having constructed some limiting cases analytically, we now examine the full parameter space by scanning through 
it, solving for G ff numerically. Our aim is to find generic trends in the behaviour of G^' and G^g , and furthermore 
restrict the parameter space by applying theoretical constraints on the behaviour of the ir field. 

As an initial illustration of the paths of gravity, i.e. the evolution of the gravitational strength G c ff (a), we consider 
the family of models that result in G Q ff.ds/GN = 1. For the uncoupled case cq — cq = our approach is to fix 
C5 = —1, and C2 by demanding 0^,0 = 0.76. There remains a set of (03,04) values that can satisfy our remaining 
condition G e g,ds/GN = 1- 

Figure [6l right panel, plots the effective gravitational strengths G^g (a) / G-® , G^\a)/G^, and G^ r ^\a) / 'Gn for 
one case giving rise to the asymptotic behaviour G e ff,ds/GAr = 1. Despite having past and future asymptotes agreeing 
with general relativity, the intermediate evolution over the redshift range z ~ (0, 10) can be strongly distinct, both 
from general relativity and among the different G G ff's. The magnitude of the features over this range is extremely 
sensitive to the values of 03,4. 

While it is instructive to analyze particular examples of model parameters, we can also learn a lot by a broader 
if shallower study. Generically we have a three (uncoupled case) or four (linearly or derivatively coupled cases) 
dimensional parameter space spanned by €3-5, plus possibly cq or cg- As before we fix C2 by requiring the present 
dark energy density = 0.76. In what follows we also fix C5 = —1 for tractability, and hence as mentioned in 
the previous subsection only have to consider the range 03,04 < 0. We systematically scan the parameter space by 
selecting 22,500 points from a uniform grid over the range (—5, 0) for each of the 03,4 terms. We discuss the Co and cq 
parameters below. For each point we evolve the background equations (|8| 9p from z = 10 6 to z = —0.99 (i.e. future 
scale factor a = 100). From the 22,500 expansion histories, we keep as viable only those that satisfy the following 
conditions for z < 500: 

(a) The energy density of the 7r field is positive definite, p n > 0. 

(b) The no-ghost condition (|P8[) is satisfied. 

(c) The Laplace inequality \E2\ is satisfied: c 2 s > 0. 

The positivity of the energy density is not strictly required at all times during the cosmological history. However, 
there are a number of reasons as to why we include it. We have already seen in ([5^1 - I 56[) that at early times the 
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FIG. 7: [Left panel] Histogram of the gravitational strength G^j /Gn for an ensemble of theoretically viable models arising 
from a uniform scan through the uncoupled Galileon parameter space, at redshift z = 1 (dark) and z — (light). At z = 1 we 
find G e ff to be clustered in the vicinity of its general relativistic limit (deviating by ~ 10%), however at late times the behaviour 
is more diverse and further from the GR limit. We observe that G e ff can be both positive and negative at late times, although 
the majority of the parameter space surveyed approached G e g > 0. [Right panel] The asymptotic value of G e g at z = —0.99 
(essentially the de Sitter asymptote), showing extremely diffuse behaviour in the future. 



no-ghost constraint is intimately tied to the p v > condition. Similarly, we must have p^ > at the present time 
since we demand that the field drives the current late time acceleration. This does not preclude the possibility that 
the energy density can be positive at high redshift, cross through p^ = at some intermediate time and then return 
to a positive value at z = 0. It is very difficult to prove that such a crossing can never occur, however we have reason 
to believe that such behavior is highly unlikely. There are multiple fixed points that the galileon field can approach 
asymptotically; in addition to the de Sitter point there are both anti de Sitter and Minkowski vacuum states, and it 
is likely that any crossing of the zero line p^ — will result in the galileon failing to approach the dc Sitter state that 
we require. Hence we impose p n > 0. 

Figure [7] plots the histograms of the results for the gravitational strengths G^g at redshifts z = 0, z = 1, and 
z = —0.99 for each of the surviving, theoretically viable runs. At z = 1 the viable models are strongly clustered in 
the vicinity of G^g /Gn > 1, with deviations from the high redshift, general relativistic (GR) value by only ~ 10%. 
However by z = 0, when cosmic acceleration is strong, gravity generically is quite distinct from GR and a diversity 
of behaviors is possible. At late times, essentially at the de Sitter state, the value of the gravitational strength is 
diverse, and takes both positive and negative values of G^ /Gn (cf. the right panel of Figure [SJ). If we were to fix 
G e s,ds/GN = 1) say, then we find much greater diversity at z ps 1 but of course uniformity in the future. 

We next perform an analogous numerical study for the linearly coupled Galileon model, expanding our parameter 
grid. We continue to fix C5 = —1, but now consider 250,000 expansion histories over the range C3 = (—5,5), 
C4 = (—5,5) and cq = (—5,5). The results are exhibited in Figure [H In this case, we no longer demand that our 
theoretical constraints are satisfied at all times, but rather only for z > 0. We do not consider z < as there is no de 
Sitter fixed point for this model, so we would have to check infinitely far to the future, and z = is a logical cutoff 
in choosing a domain of applicability. 

The distribution of G^ is similar to the uncoupled case - recall from Sec. HVBI that generally the cq contribution 
has to be subdominant at early times. However the later time behaviour is more diffuse. This is due to the fact that 
we now permit G c s to diverge for z < 0. Such divergences can occur for any z < 0, and the value of G e g at z = is 
particularly sensitive to the existence and redshift of these features. 

Finally, we repeat our analysis for the derivatively coupled Galileon model, expanding our parameter grid to include 
cq 0. We evenly distribute 250,000 points amongst the three dimensional space C3 = (—5,5), C4 = (—5,5), cq = 
(—5,5). The resulting histograms of theoretically viable models are exhibited in Figure [5] We observe qualitatively 
similar behaviour to the uncoupled model at both z = 1 and z = 0, although we find considerably more diffuse 
behaviour at both z = and z = 1, with a higher propensity towards negative asymptotes. 

Over the course of our numerical study, we considered the range cq = (—5, 5). However we found no models which 
were deemed viable for cq < 0. The theoretical constraint that causes the derivatively coupled model to fail in this 
regime is the Laplace condition of Appendix [Ej > 0. If we evolve our equations back to radiation domination in the 
early universe, whenever the cq term is dominant then = — 1/3 + 0(fl 7T (z)), indicating a classical instability [56| . 
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FIG. 8: Histogram of the gravitational strength G^g /Gn for an ensemble of theoretically viable models arising from a uniform 
scan through the linearly coupled Galileon parameter space, at redshift z = 1 (dark) and z = (light). We do not here consider 
models that diverge in the future to be unviable, so greater diversity is exhibited compared to the uncoupled model. 





FIG. 9: Histogram of the gravitational strength G^ /Gn for an ensemble of theoretically viable models arising from a uniform 
scan through the derivatively coupled Galileon parameter space, at redshift z = 1 (dark) and z = (light) in the left panel, 
and z = —0.99 in the right panel. We observe similar, although more diffuse, behaviour in the model with cg / compared to 
the uncoupled case. 



Although we do not consider cosmological perturbations at such high redshifts in this work, the Laplace condition 
is applicable in the early universe since to leading order the metric perturbations decouple then, and we can simply 
consider perturbations of the Galileon field evolving on the background. In this case, c 2 = —1/3, but even during 
matter domination, whenever /9 ra( j > we expect, and observe, this instability to be present, leading to the absence 
of viable models in our scans that differ substantially from the uncoupled case. An important consequence is that 
since the cq term is always of order H 2 larger than the C2 term, a model with only C2 and cq, such as that of 44 1 
(even generalized to arbitrary C2 so as to enable a de Sitter future state), is ruled out. Also, the fact that only cq > 
models survive our constraints implies that this term cannot dominate p n during matter and radiation domination. 

In Figure[in]we exhibit the evolution of G^g(a) for a "typical" viable model - that is, a model near the peak of the 
histograms for z = 1 - for each of the uncoupled, linearly coupled, and derivatively coupled cases. We note a common 
deviation from general relativity at z w 10, a reapproach to GR, and then the strong departure as cosmic acceleration 

occurs. Evolving to the future beyond z = 0, G^g approaches a positive, constant value in the de Sitter asymptote for 

the derivative and uncoupled cases, while in the linear case will continue to grow due to the explicit dependence 
on y. 

The nonmonotonic behaviour of G^g at z ~ 10 is a common feature over the parameter space considered, and 
its origin can be traced to an interaction between the different c n terms. Although the C5 term, say, dominates at 
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FIG. 10: Typical examples of [Left panel] G^(a) and [Right panel] w(a) for viable models in the uncoupled (solid red), linearly 
coupled (long dash black), and derivatively coupled (short dash blue) classes. For each of the cases we fix C3 = —5.73, C4 = —1.2 
and C5 = —1, and adjust the parameter Co or cg such that the model exists in the vicinity of the peak of that case's G e s(z = 1) 
histogram (specifically Co = 0.05 or cg = 0.4). The behaviour of uncoupled and derivatively coupled models is qualitatively 
similar, with both gJ^ functions exhibiting a late time approach to a constant positive asymptote. 

early times, the 0^4 contributions grow relative to this and at some redshift they become comparable. This leads to 

cancellations in the denominator of GjL and account for the observed bump. The height and redshift of this feature 
depend on the values of C2-5 and also p v {zi); for example for lower initial 7r density, the hierarchy between the C3_5 
terms is less pronounced, leading to a feature at earlier times. Gcnerically, then, the paths of gravity in Galileon 
gravity are quite intricate and rich. 



VI. CONCLUSIONS 



Galileon gravity is a well defined theory with a rich phenomenology. We explored the predictions for expansion 
history giving rise to current acceleration, despite the absence of any cosmological constant or potential, and found 
early time scaling and tracing solutions and late time de Sitter states. Equally important is the behavior of the 
inhomogencities in the Galileon field, and we discussed requirements for a sound theory in terms of the no-ghost 
condition, Laplace stability, and positivity of energy density. By creating a large ensemble of expansion histories, we 
were able to use these theoretical conditions to efficiently rule out large regions of the multi dimensional parameter 
space. In this work we focussed solely on the dynamics of scalar perturbations; if we include tensor modes then we 
expect even more stringent bounds on the Galileon parameters (see for example [53j). Confronting the model with 
data to further constrain the parameters will be the subject of future work. 

The modifications to gravity induced by the Galileon terms can be characterized by gravitational strengths (effective 

Newton's constants) entering the Poisson equations for the various combinations of metric potentials, G^g , G^\ and 
G (0+V-)_ Wc der i V ed expressions for these and solutions in the early time and late time limits. Considering special 
cases and general scans over the parameter space, we presented the "paths of gravity" , the evolution of the effective 
gravitational strengths over cosmic history. 

Galileon models do not have simple, monotonic evolution of gravity. One special feature discovered is that in the 
de Sitter limit the "matter gravity" G^g entering into growth of structure and the "light gravity" G^ + ^ entering 
into light deflection are identical (unlike in many other theories), although generally different from general relativity. 
Gravity can suddenly become stronger at some epoch in 1 < z < 1000, and then be restored to general relativity; it is 
intriguing to consider the possible effects on structure formation such as enabling early formation of massive clusters. 
The observed phenomenology of G^g at intermediate redshifts will significantly increase the (linear) power of the 
matter perturbations at late times, and this will be used to severely constrain the uncoupled Galileon in an upcoming 
study. 

Models which include a coupling between the Galileon field and the energy-momentum tensor, or cquivalcntly the 
field and the Ricci scalar in the case of linear coupling or the field kinetic term and the Einstein tensor in the case 
of derivative coupling, are severely constrained by conditions for soundness of the theory. In particular, we find that 
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each coupling term must be subdominant to uncoupled Galilcon terms in the early universe, and generally small. 
Deviations from the uncoupled case therefore only tend to show up at late times. For example, the linearly coupled 
case does not have a true de Sitter fixed point. 

Galileon gravity offers a rich variety in both the expansion behavior (e.g. equations of state that can track the 
matter component and give rise to early dark energy) and gravity behavior (spikes of deviations from, and restoration 
to, general relativity). It provides a theoretically viable alternative to Einstein gravity, motivated by higher dimension 
geometric theories and protected by symmetries, without many of the naturalness issues of the cosmological constant 
or scalar field potentials. Galileon cosmology's implications for observations at high and low redshift make it an 
excellent theory to explore, and to test as an origin for cosmic acceleration and an extension to Einstein gravity. 
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Appendix A: Jordan Frame 



In the main body of the paper we choose to work in a frame in which there is no direct coupling between the it field 
and standard matter, instead the coupling is exhibited directly to the metric quantities. To accomplish this we here 
consider the transform between the Einstein and Jordan frames in the weak field limit (see also [42|). We consider 
two couplings throughout this work, and we treat them separately here. Beginning with the linear coupling, we write 
9[iv — f)u,v + hn V and the Lagrangian to quadratic order in h^ v is given by 



C 



M 2 pl 



h» v 6G 



1 5 1 

i — 1 r 



(Al) 



where SG^ is the linearized Einstein tensor and T — rf^T^^ where is the energy momentum tensor of matter. 
In this linearized limit we can perform the transformation /i„ y — > hn V = hu, v + Icqit jM p \r\u. v to remove the Cq coupling 
term. Under this transformation, neglecting any boundary terms, the linearized Einstein contribution in Eq. (IA1I) . 
the first term on the right hand side, introduces three additional terms, two of which can be absorbed into the C2.3 
parameters. In terms of h^ v and 7r, the Lagrangian can be written as 



£ = 



M 2 pl 



1 5 1 - 

h^SG^ + 2c M pimi ^8G^ + - ]T c t C t + -hTT^ , 



(A2) 



where Ci are the redefined constants, and 5G^ V is the linearized Einstein tensor involving on h^. Having constructed 
the linearized Jordan frame limit, we promote the action to its full, nonlinear form: 



S 



-gd x 



l-2c 



M p i 



_C3_ 

A/ 3 



(d7r) 2 D7r 



C4C4 c 5 £ 5 



(A3) 



Similarly, we can repeat the above procedure with the derivative coupling cqT^ d^Trd v Tr. Specifically, we write the 
linearized Lagrangian as 



C 



2 ^ 

i=l 



CiCi 



M pl AP 



(A4) 



Performing the transformation /i M „ — > = + 2{cQ/M p \M 5 )d tl ird v iT, we remove the explicit coupling to matter 
at the expense of introducing additional mixing terms between h^ v and w. The transform introduces a term that can 
be absorbed into C4, and a new term of the form d^nd v n&G ilv . We find 



M l\ - - Moi - 1 1 - 

C = -^h^SG^ + 2c G j^d^7rd v 7r5G flu + - £ cjA + -hTT^ 

i=i 



(A5) 
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Promoting to a full, nonlinear action, we arrive at 



S = J d 4 x y/—g 



M pl R C 2 2 C 3 



2 W - j^^fU-K 



CjC A _ c 5 £ 5 _ M p i 
2 2 M 3 



(A6) 



In regions of high curvature, the equivalence between the Jordan and Einstein frames considered here will no longer 
be valid, however we work with the Jordan frame quantities throughout. 

Appendix B: Field Equations 

To obtain the Einstein and scalar field equations we vary the action of Eq. © with respect to the metric g^ v and 
Galileon held 7r, obtaining 



M 3 



C-2 



2 

M pl M 3 c 



R + [-4(D7r) 3 - 8V Q V t/ 7rV l/ V A 7rV A V Q 7 r + 12(D7r)V Q V (3 7rV Q V /3 7r + 2D7rV Q 7rV Q 7ri? 



4M 3 



-4V a 7rV^7rV a V /, 7r J R + 8(D7r)i?^V /1 7rV I ,7r - 4V a 7rV Q 7rV A V T 7ri? A(T - 16V a 7rV a V <T 7ri? ap V p 7r 
-8V Q 7rV /3 7rV (3 V (T 7r J R Q '" 3<T ] 



pi 



where the C2-5 contributions to the Galileon energy momentum tensor, T^ 2 5 \ were first calculated in We repeat 
them here for completeness, and add the cq contribution. 
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7 1 ( C2 ) 



rp(c 4 ) 



T'Ccb) 



t( c g) 



C2 



[2V M 7rV v 7rD7r + .g^V^V"^) 2 - 2V (/1 7rV„) (dir) 2 ] 



_£3_ 

~~ M 3 

= [4D7rV a 7rV„7rV Q V z ,7r + 4D7rV a 7rV l/ 7rV Q Vu7r - 2D7rD7rVu7rV l/ 7r + 2D7rV ct 7rV Q 7rV A1 V^7r 



(Bl) 
(B2) 



j^j-q L — '" T t*" T M ^ 1 — a ^ M — ^^^^ — M ^ 1 — a M 

+4V Q 7rV A 7rV Q VA7rV Al V I ,7r - 4V"7rV A 7rV Q V AI 7rV A V^7r + 2V A1 7rV„7rV Q! V /3 7rV Q V^tt 

-2V Q 7rV Q 7rV AI V, 3 7rV I ,V /3 7r - 4V Q 7rV t ,7rV Q V (T 7rV ff V Al 7r - 4V Q 7rV AI 7rV Q V CT 7rV CT V^7r - g /lt/ (n 7 r) 2 V a 7rV a! 7r 
-4.9 Mt/ n7rV Q 7rV A 7rV Q V a tt + 45^V Q 7rV LJ 7rV Q V A 7rV A V"7r + p^V Q 7rV Q 7rV A V (T 7rV A V' T 7r 



-V M 7rV„7rV a 7rV a 7r.R ■ 



9/tV 



V„7rV a 7rV>7rV A 7ri? - 2V„7rV Q 7ri? AM V !/ 7rV A 7r - 2V„7rV Q 7ri? A „V„7rV A 7r 



l Ai ,V p 7 



_£5_ 

M 9 



-i? Miy V Q 7rV Q ^V A 7rV A ^ + 2< ?/lt/ iZ CTW V< T 7rV a, 7rV C( 7TV Q 7T - 2V a 7rV a 7n ?K i^ tn , w V< T 7rV" ; 7r 
(□tt) 3 V^ttV^tt + (Dtt) 3 V Q 7rV Q 7r 5A1!/ - 3 (Dtt) 2 V a 7rV a 7rV (U V i/ 7r 



(B3) 



-3 (Dtt) 2 V p tt (VV^ttV^tt + VV^ttV^tt) + 3 (Dtt) 2 V Q 7rV (3 7rV Q V /3 7r.g AIl/ 

+6 (Dtt) V Q 7rV Q 7rV M V (T 7rV' T V I ,7r - 6(D7r)V Q 7rV /3 7rV Q V, 3 7rV M V 1 ,7r - 3(D7r)V Q V /3 7rV ct V^ttV^ttV^tt 

+6(D7r)V a 7rV /3 7rV Q V Al 7rV, 3 V ly 7r + 6(D7r)V Q 7rV Q V,37r (V^V^ttV^tt + V^V^ttV^tt) 

3 

-3(D7r)V C[ 7rV Q 7rVAV (T 7rV A V ,T 7r.g^ - 6(D7r)V Q 7rV /3 7rV Q V „-kV pVitg^ - -{Utt)V a nV a irV y iiR 
+3(D7r)V Q 7rV Q 7rV /3 7r {R^V^ + Rp v V ^) - 3(D7r)V Q 7rV Q 7rV A 7rV /3 7ri? A , 35Ait/ 

+3(D7r)V Q 7rV Q 7rV A 7rV ,T 7ri?^ A + 3V Q 7rV Q 7rV A V ,T 7rV A V a 7rV Al V I ,7r - 6V Q 7rV Q 7rV A1 V ,T 7rV (T V A 7rV A V„7r 
+6V Q 7rV /3 7rV Q V / 37rV A V M 7rV A V ly 7r + 6V Q 7rV' 3 7rV Q V CT 7rV CT V / 37rV A1 V t/ 7r + 2V (9 V' T 7rV (T V A 7rV A V p 7rV M 7rV t ,7r 
+3V Q V^7rV Q V / 37rV A 7r (V A V M 7rV„7r + V A V„7rV A1 7r) - 6V p 7rV p V' T 7rV (7 V A 7r (V A V M 7rV„7r + V A V„7rV AI 7r) 
-6V7rV p V A 7rV (T (V A V M 7rV a V I ,7r + V A V t/ 7rV' T V M 7r) + 2V Q 7rV Q 7rV CT V A 7rV A V K 7rV K V ,T 7r ftu , 
-3V Q 7rV /3 7rV Q V /3 7rV A V' T 7rV A V ff 7rc/^ + 6V p 7rV K 7rV p V (T 7rV' T V A 7rV A V K 7r ffAII , 



3 3 

--V a 7rV a 7rV CT 7r (VctV^ttV^tt + V^V^ttV^tt) R - -V a 7rV Q! 7rV A 7rV' T 7rV A V CT 7riJ5( 



+3V a 7rV a 7rV A 7rV <T 7rV A V [T 7r J R Atly + 3V a 7rV Q 7rV A 7rV <T 7ri? ACT V M V i ,7r + 3V a 7rV a 7rV A V IT 7ri? A(T V p 7rV ! ,7r 
-3V Q ^V Q 7rVVV A V ff 7r {R^VvK + i? CT ,V^) - 3V Q 7rV Q 7rVVV CT 7r (R Xlt VuV a ir + i? A „V M V CT 7r) 
~3V a TrV a TTV 7 R aX (V A V M 7rV„7r + V^ttV^tt) + 6V Q 7rV Q 7rV <7 7rV CT V A 7ri? AK V K 7r ffM1/ 
-3V Q ^V Q 7rV ff 7rV A V K 7r {R^Vyir + i?„ ACTK V M 7r) + 3V Q 7rV Q 7rV ff 7rV A 7r (R^ aXK \7 K \7 „ir + A^VV^tt) 
-3V Q ^V Q 7rV ff 7rV a VVV K 7r (R„ XvK + R uX „ K ) + 3V a i:V a nV a nV\V K V T nR aKXT g^] (B4) 



M p i 
A/ 3 



c G [^DTi-Ck - 2D7rV A1 V„7r + 2V M V A 7rV„V A 7r - 5 A1! ,V A V Q 7rV A V Q 7 r] 



(B5) 



M 3 



CG 



+2^|c G [i? A „V A 7rV M 7r + i? v V A 7rV„7r - ^i?p A V p 7rV A 7r + iT^WV^Tt] 



Appendix C: Perturbed Field Equations 



To obtain the perturbed field equations in the sub horizon limit, taking the Newtonian gauge and defining Sy = 
6n /M p i, we linearize the perturbed i) and (0,0) Einstein, 7r dynamical, and fluid conservation equations in 5y, 
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the metric potentials <f> and ip, and the density perturbation S m . We obtain the following equations 
- [1 - 2c y] (V 2 V> - V 2 0) = -2c V 2 Sy 



-C4 



1 - \ iJ 4 x 4 - 

6H 3 x 2 ( #' a; + Hx' + -.Hz J V 2 <5y + — — V 2 c, 



-H 4 x 4 V 2 i> 



+c G [2 (HH'x + H 2 x' + H 2 x) V 2 5y + H 2 x 2 V 2 (f) + H 2 x 2 V 2 ip] 

+c 5 [3H 6 x 5 V 2 ^ - 3H 5 x 4 (Hx' + H'x) V 2 <p + H 5 x 3 (12 Hx' + 15H'x + 3Hx) V 2 6y] (CI) 



\ [1 - 2c y] V 2 <^ = 



S ni + 2c 3 H 



2c 



"CG 



H 2 M 2 r 

2 

AH 2 ^V 2 5y-2H 2 %V 2 c, 
a 1 



a- 



C5 



or 

15i7 6 x 4 



- 2 _ (\2H 4 x 3 



\7 2 8y- 



3H 4 x 



4^4 



2x - V 2 <; 



(C2) 



c 3 HH'x + c 3 H 2 x' + 2c 3 H 2 x 



C2 



V 2 5y - c 3 H 2 x 2 V 2 ip = -c [V 2 4> - 2V 2 0] + c 4 [-6H 4 x 3 V 2 ip 

+ (6H 3 x 2 (Hx' + xH') + 2H 4 x 3 ) V 2 <j) - (\2H 3 x (Hx' + xH') + 13ff 4 T 2 + 6H 3 x 2 H') V 2 5y] 
+c G (2HH' + ZH 2 ) V 2 Sy + 2c G H 2 xV 2 tp - 2c G (S (Hx' + H'x) + H 2 x) V 2 

(l8H 6 x 2 x' + 30H 5 x 3 H' + 12H 6 x 3 ) \7 2 Sy + ^-H Q x 4 V 2 ^ - (3H 6 x 4 + l5H 5 x 4 H' + 12H 6 x 3 x') V 2 



-C5 



H 2 5'^ + HH'S' + 2H 2 S' = -^V 2 ^ 



(C3) 
(C4) 



where V = \7/Hq, p m is the matter density, and <5 rn = 5p m / p m . We have assumed that the quasistatic approximation 
holds at all times: 5y,4>,4> <C V 2 5y, V 2 0, V 2 ?/>. Note that in Eq. (|C1[) the "slip" between and ip is not sourced by 
c-i or C3. The above equations are in agreement with those derived in [30| in the uncoupled case cq = cq = 0. 

Equations (|C1[) and (|C3[) can be used to eliminate V 2 ch/ and V 2 ^ from the Poisson equation (|C2[) . Hence we 
can describe the effect of the Galileon on subhorizon density perturbations through two functions of the background 
quantities H,x,y. We define these functions G { £ and G^ in Eqs. (22) and ([25]), and = + G ( ^>]/2. 



Appendix D: No-Ghost Condition in Detail 



For a given coupling, the action © describes a model with five free parameters: the magnitude of the coupling 
strength and the kinetic terms (02, C3, C4, C5). However, one must be careful to restrict our analysis to the parameter 
space in which the model is theoretically viable. Of particular importance is the absence of propagating ghost degrees 
of freedom, that is degrees of freedom whose kinetic contribution to the Hamiltonian is negative, making it unbounded 
from below. 

In Q a "no-ghost" condition was derived for models of the form 



d 4 Xy 



\f{k)R + i/ a (7r, X) + C(tt)(cV) 2 & 



(Dl) 



where X = — 8^8^/2, and F(ir), f 2 (TT,X) and ((tt) are arbitrary functions. The no-ghost condition corresponds 
to @ 



(24C1/7T - 8C> 2 + h,x + f2,xxK 2 ) Fn 2 + 3 (P - 2^^ > 



(D2) 



The action (j-Pip is equivalent to ours with C4 = C5 = cq = cq = 0; in this case we can use the above condition as 
a constraint on our parameter space. When we introduce the couplings co,g and more general kinetic terms C4 ; 5 the 
analysis must be redone as follows. 

To construct the no-ghost condition, we utilise four equations: the tt equation of motion and the (0,0), trace, and 
(i,j 7^ i) Einstein equations. To deduce whether any propagating degrees of freedom exhibit ghost behaviour, it 
suffices to look for the sign of the coefficient of the cj> and by terms. Our approach is to note that the 7^ i) 
equation gives an algebraic relation between Sy, tp and and the (0, 0) Einstein equation contains only V 2, 0, V 2 
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and ~S7 2 5y terms (and first order time derivatives, which are not important in what follows). Therefore combining 
these equations yields a relationship between V 2 and ~S7 2 8y. The ir equation of motion and trace Einstein equation 
contain the terms Sy, (j>, V 2 f/>, W 2 (j> and V 2 <5y, and therefore by eliminating V 2 ip and using our relationship between 
V 2 <5y and V 2 we can construct equations of the form 



Ax + BV 2 x = C . (D3) 

Here A, B are diagonal 2x2 matrices that are functions of background quantities such as tt, H etc. and x is the two 
dimensional vector containing Sy and (f>. The matrix C contains lower order time and spatial derivatives of x and is 
unimportant in calculating the no-ghost condition. For the simple case of a minimally coupled, canonical scalar field 
containing the c 2 term only, the no- ghost condition corresponds to C2 > 0, and therefore we simply need to ensure 
that A preserves the same sign conventions. 

We begin with a simple case taking 02,3 only. The 7^ i) equation gives ip = </>, and the (0,0) component of the 
Einstein equations yields the relationship 



^=^mM^V + ... (D4) 

which will not be necessary in what follows. The . . . denote lower derivative contributions to the equations that will 
not be relevant to the stability condition derived here. 
The trace and 7r equations are 

(60 - 2V 2 0) = (67r 2 Sy - 2rr 2 V 2 Sy) + . . . (D5) 

-Sy + V 2 5y) + 6c 3 HnSy - 2c 3 (tt + 2Hn) V 2 5y - c 3 n 2 Uij) - V 2 <p] + • • • = (D6) 



2 

Eliminating <p, we find that the coefficient of the Sy term is given by 



which must be less than zero to ensure that perturbations of the tt field have the correct sign. This condition is in 
agreement with Eq. (|_P2[) . 

We now expand our approach to include the 04,5^ terms in the action. Performing exactly the same steps as above, 
we find the following inequality for the no-ghost condition 

c 2 -o -0 -49 -« , (3c 3 H 2 x 2 + 6c G H 2 x- 18c 4 i? 4 .T 3 + f c 5 H 6 x A -3c ) 2 , . 

This must be satisfied at all times during the cosmological evolution to ensure that at the level of linear perturbations 
around the cosmological background Sir has a kinetic term that contributes positively to the Hamiltonian. Note that 
the no-ghost condition (|Z38[) can be written in terms of the Ki functions as 



«3 + |^<0. (D9) 

Returning to our solutions for the linear coupling case, Eqs. (|45l) - (|47p valid during matter domination, we can show 
that for parameters that give a positive energy density > 0, the condition (|_D8j) is not satisfied. To see this, we 
note that at early times the no-ghost condition is well approximated by the first term K2, 

- y + 6c 3 H 2 x + 3c G H 2 - 27c4,H 4 x 2 + 30c 5 H 6 x 3 < . (D10) 

Consider the various dominating pairs, starting with C5, cq- To ensure that > 0, we are forced to take cq > as 
Aq must be real. This in turn forces C5 > 0, which violates Eq. (|£>10j) . We find the same dilemma for the C4 solution 
and so on until we reach c-i . Note that we have not shown that a ghost is generically present for the linearly coupled 
model whenever p^ > 0, but rather just for the particular solutions found in llV B[ In fact during radiation domination 
the no-ghost condition is generically satisfied, as in the uncoupled case. 

The no-ghost condition derived here is applicable to linear perturbations on an FRW background. A more com- 
plete analysis should be undertaken in regimes where the non-linear nature of the derivative self couplings becomes 
significant. 
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Appendix E: Laplace Instability 

In the previous section we considered the sign of the Sy term in the linearized perturbation equations, and how 
the no-ghost condition can be used to place constraints on the Galileon parameters. In this section we consider the 
existence of a second instability that must also be avoided - the Laplace instability [2!| (see also [H?} for a very 
complete discussion of instabilities). This is a condition on the sound speed of the it field perturbation, obtained by 
linearizing the Einstein and tt equations and eliminating the metric potentials to obtain a wave-like equation for Sir. 
Specifically, we combine the (0,0), trace, and ^ i) Einstein equations with the perturbed tt equation of motion, 
keeping all terms containing second order time and spatial derivatives of the fields. We find the following expression 
for Sir: 

3 k§\ / 2k3«| + 2k 2 kq — 4K1K4K5 

2 K4 J V 2k 4 



"a + i^ **+ — T. — V*frr = ... (El) 



Again, . . . represents terms unimportant to the stability argument. The negative defmitcness of the Sn coefficient 
is the no-ghost condition, and the negative defmitcness of the ratio of the V 2 8ir and Sir coefficients is the Laplacian 
stability condition. This ratio is the negative of the sound speed squared, 



_ 4kik 4 k 5 - 2k 3 kI - 2k|k 6 

c s = 77 1 o 2 \ ■ y^ 1 ' 

K4 (ZK4K2 + JK5) 



The Laplace stability condition is c 2 > 0. If this condition is violated, imaginary frequency behaviour would follow, 
leading to exponential growth of the tt field perturbation. (It is conceivable that such behaviour is benign if the 
timescale associated with the growth is suitably large, however we opt for a conservative approach and enforce c 2 > 
at all times during the cosmological evolution.) 

At early times when n n <C 1, the sound speed simplifies to c 2 = —k§/k,2- At late times, for the uncoupled case we 
can derive a relatively straightforward asymptotic de Sitter form for the sound speed, corresponding to 

(14 + 3C - AD) (416 + 108C + 9C 2 - 32D - IQD 2 ) 

54 (8 + 3C - AD) (20 + 3C - AD) ~ ^ ' 

where C and D are defined in Eqs. (jrJ2"j) . (|63l) . 

In the main body of the text, we impose the positivity of p w and the no-ghost and Laplace conditions for all z < 500. 
There is an additional condition that one might impose, that the Sir field must avoid "superluminal" behaviour c 2 > 1. 
However, if such a condition were violated it is not necessarily the case that the model is ruled out; see for example [58[ 
for a detailed discussion of superluminal propagation in a different class of scalar fields. In addition, a more complete 
analysis should take into account the full perturbation equations when constructing the scalar field dispersion relation, 
including effective mass terms of the form H 2 Sy. For these reasons we note the potential existence of this constraint 
but do not impose it when scanning over the parameter space for viable models in Sec. IV Dl 
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